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Advancing the debate on the possible advantages of quantum effects in biomolecular systems requires illustra- 
tions of biomolecular prototypes that can exhibit and use non-classical phenomena to aid performance. Here 
we investigate non-classical fluctuations in molecular dimers present in photosynthetic light-harvesting anten- 
nae. We show that exciton-vibrational resonances enable non-classical occupation fluctuations of the driven 
collective high-energy vibration at room temperature. This phenomenon is a manifestation of strong mod- 
ulation of low-phonon occupation levels through the transient formation of hybrid exciton-vibration states. 
For the prototype considered, direct quantitative relations seem to hold between non-classicality and optimal 
exciton population transfer. 



The experimental investigation of quantum effects in 
photosynthetic energy transfei^^ has inspired exciting 
research towards asserting the relevance of non-classical 
behaviour in the biological function of the molecular com- 
ponents of living organisms. Theoretical works have 
made emphasis on the quantumness exhibited by elec- 
tronic degrees of freedonP and on the key role of envi- 
ronmental fluctuations^^. In order to understand the full 
range of possible non-classical phenomena present dur- 
ing room-temperature dynamics of electronic excitations 
in light-harvesting antennae it is important to expand 
the analysis to the molecular motions that play a fun- 
damental role in driving exciton transport^. Res earch 
in this direction is timely given that experiment J 10 ! 11 -! 
and theoretical worka^SHISl are highlighting the impor- 
tance of concerted electronic and vibrational dynam- 
ics in energy transfer. Moreover, recent experimen- 
tal techniques able to manipulate vibrational states^^oi 
and probe their quantum properties may indeed pro- 
vide experimental access to the non-classicality of vi- 
brational modes in biomolecular systems. Here we con- 
sider a prototype molecular dimer with an structured 
exciton-phonon interaction and investigate non-classical 
behaviour of molecular vibrations that steer energy dy- 
namics and exciton coherences in this molecular unit. 
While previous studies have considered coherent exciton- 
vibration dynamics in molecular dimers^ 4 * 17 * 21 !, it is un- 
clear how this deviates from classical (wave-like) coherent 
behaviour and if so whether there are quantitative re- 
lationships between non-classicality and effective energy 
transfer. 

The prototype of interest consists of two chromophores 
coupled by weak electronic interactions and subject to 
the influence of both a low-energy thermal background 
and a discrete high-frequency intramolecular vibration. 
The key aspect of this molecular vibration is that its 
energy matches closely the exciton gap but it is much 
larger than the thermal energy scale. Dimers of this type 
naturally occur in several light-harvesting protein*^!] 



a ) Electronic mail: a.olaya@ucl.ac.uk 



and are of particular interest as they represent sys- 
tems where, although counter-intuitive, exciton coher- 
ence beating has been probecP^. Insights into the im- 
portance of resonances between exciton energy mismatch 
and vibrational frequencies for efficient energy trans- 
fer can be gained from Forster theory^. However, the 
wider implications of such reso nance s for optimal spatio- 
temporal distrib ution of energy^ 2 ^, for manipulation of 
exciton coherenc e! 13 * 14 * 16 * and for non-adiabatic excitation 
dynamical have only started to be clarified. An is- 
sue that remains to be tackled, is the extent to which 
such resonances can expose the non-classical features 
of the molecular motions activating transport. In the 
dimer investigated here, commensurate energies of ex- 
citon gap and discrete vibrations imply that after in- 
put of exciton population, effective energy distribution 
proceeds by driving the collective relative displacement 
mode out of thermal equilibrium towards a non-classical 
state. Non-classicality of the state of this relative motion 
is unambiguously indicated by sub-Poissonian occupation 
statistics 2 ^ or alternatively by strong modulation of low- 
phonon occupations 2 and is a manifestation of transient 
collective coherent behaviour between excitons and vibra- 
tions. Low-energy thermal fluctuations stabilize energy 
distribution to the lowest exciton state before destroy- 
ing non-classicality of the collective mode assisting the 
process. Meanwhile, coherent interaction between elec- 
tronic degrees of freedom and such a mode manifests it- 
self in the beating of excitonic coherences with frequency 
components that vary depending on the coupling to the 
thermal bath. When time-averaging over a relevant time 
scale, we find direct quantitative relations between pop- 
ulation transfer and the degree of non-classicality as in- 
dicated by the average of negative values of the Man- 
del Q— parameter^. This relationship holds even when 
the electronic system is initialized in a statistical mixture 
of exciton states. Hence, efficient photosynthetic energy 
transport could take advantage of non-classicality. 

The field of quantum optics has developed a solid 
framework to quantify the non-classicality of bosonic 
fields^. It therefore provides useful conceptual and 
quantitative tools to investigate non-classicality of the 
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harmonic vibrational degrees of freedom of interest in 
this work. From the perspective of quantum optics, non- 
classicality arises if the state of the system of interest 
cannot be expressed as a statistical mixture of coher- 
ent states defining a valid probability measure, this then 
leads to non-positive or highly singular values of the 
Glauber-Sudarshan P-function or the Wigner function 29 . 
However, experimental characterization of these distri- 
butions is not always a straightforward task. Alterna- 
tively, it is of interest to investigate signatures of non- 
classicality. For a single-mode field, a sub-Poissonian oc- 
cupation number distribution quantified by the Mandcl's 
Q— parameteJ^fil is a signature of deviations from classical 
behaviour. The Mandcl Q— parameter characterizes the 
departure of the occupation number distribution P{n) 
from Poissonian statistics through the inequality 



where (n) and (n 2 ) denote the first and second moments 
of the bosonic number operator n respectively. Van- 
ishing Q indicates Poissonian number statistics where 
the mean of n equals its variance (such is the case of 
a coherent state of light). For a chaotic thermal state 
one finds that Q = (n) > indicating that particles 
are 'bunched'. The departure from classical (wave-like) 
behaviour is characterised by Q < or particle num- 
ber distributions where the variance is less than the 
mean, as is the case for Fock states. The occupation 
fluctuations associated with a sub-Poissonian distribu- 
tion originate from the quantized nature of the field and 
cannot be described in terms of an underlying classi- 
cal stochastic process. A criteria that directly connects 
non-classical occupation fluctuations and the negativity 
of a quasiprobability distribution is an inequality indi- 
cating modulation of nearest number occupations^. In 
this work we use the above framework to quantify non- 
classicality of the undamped vibrational modes support- 
ing exciton transport in the molecular dimer of inter- 
est. We consider each chromophore with an excited elec- 
tronic state of energy strongly coupled to a quantized 
vibrational mode of high-frequency (w v ib > KgT) de- 
scribed by the bare Hamiltonians H c \ = 2 £,cr^"a~ 

and -ff V ib = uj V ib{b\bi + b\b2) respectively. Inter- 
chromophorc coupling is generated by dipole-dipole in- 
teractions of the form f/d-d = V(vt a 2 + °2 °T)' The 
electronic excited states interact with their local vi- 
brational environments with strength g, linearly dis- 
placing the corresponding mode coordinate, Hel-vib — 
d12i=i 2 ^t^TW + I n tne a bove b\(bi) creates (an- 
nihilates) a phonon of the vibrational mode of chro- 
mophore i while af creates or destroys an electronic 
excitation at site i. The eigenstates \X) and \Y) of 
H c \ + -f/d-d denot e exciton states with energy splitting 
given by AE = J (e\ — E2) 2 + W 1 . Transformation into 
collective mode coordinates shows that only the mode 
corresponding to the relative displacement of the modes 



with the bosonic operators feW = (pi — b 2 )/\/2, cou- 
ples to the excitonic system. The centre of mass mode 
decouples from the electronic degrees of freedom. The oc- 
cupation of the relative displacement mode of frequency 
w v ib is denoted by n and the non-classical properties of 
such mode is the subject of interest in this article. 

Our prototype system corresponds to the central dimer 
of the phycoerythrin 545 (PE545) antennae found in 
cryptophyte marine algae^l. Exciton coherence beating 
involving this dimer has been probecP and a structured 
spectral density has been phenomenologically determined 
via linear spectroscopy^. The electronic parameters for 
the dimer taken from reference^ are E\ — 19574cm -1 , 
£2 = 18532cm" 1 and V = 92.2cm -1 . The intramolecu- 
lar mode that is separated from the set of vibrational 
degrees of freedom has frequency u vih = 1111cm -1 
being quasi-resonant with the exciton energy splitting 
AE = 1058.2cm -1 . The coupling strength of each 
chromophore to the this mod e is g — WvibVO.0578 = 
267.1cm -1 obtained froir i 3 ° l 31 J, and the thermal energy 
scale at room temperature is /cbT ~ 200cm -1 . We con- 
sider the influence of the low energy thermal fluctua- 
tions via a thermal bath Hg described by a continuous 
distribution of harmonic modes, each linearly and inde- 
pendently coupled to the excited electronic states H e \ 
with coupling strength characterised by a Drude spectral 
density J{oS) — (\VL c uj /ir) / (uj 2 + fi c 2 ) 2 with reorganiza- 
tion energy A = 110cm -1 and cut-off frequency Q c — 
lOOcm -1 ^, hence the strength of the system-bath inter- 
action is parametrized by \J Since J(w V ib) <C g 2 , 
this does not over-account for the influence of the res- 
onant mode on the electronic system. Importantly, al- 
though the electronic coupling V is the smallest energy 
scale in the prototype dimer, the exciton dynamics and 
phonon fluctuations here investigated will not be prop- 
erly described by Forster theory. In order to accurately 
account for the effects of the thermal bath we have ex- 
tended the hierarchical expansion of the dynamics^^H 
to explicitly include the quantum interactions between 
electronic excitations and the selected mode. This ap- 
proach allows us to account accurately for the effects of 
the low-energy thermal bath. In calculating dynamics we 
consider an initial state, in which the upper energy exci- 
ton eigenstate is rapidly excited while both high-energy 
quantized vibrations and low-energy thermal modes are 
initially in thermal equilibrium at room temperature. 
Truncating the quantized modes at an occupation Fock 
level M = 6 adequately describes both reduced dynamics 
of the collective quantized mode and the electronic dy- 
namics. A description of the numerical methodology can 
be found in the appendix. 

Hamiltonian dynamics and non-classicality.- We first 
illustrate the quantum coherent dynamics of excita- 
tion when coupled to the undamped quasiresonant high- 
energy vibrational modes in the absence of the low-energy 
thermal bath. This allows us to identify relevant time 
scales in the evolution of a state with no initial exci- 
tonic superpositions: p(to) — \X) (X\ ® g^ h which writ- 
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FIG. 1. a, Representation of the pigments and protein envi- 
ronment of the PE545 complex (Protein Data Bank ID code 
lXGO^). The central dimer is highlighted. Also shown are the 
energy levels of the exciton- vibration states used to describe 
energy transfer \X,n). b, Coherent dynamics of the exciton- 
vibration dimer. Oscillatory population of the lowest exciton 
population pYy(t) (thick blue curve) is accompanied by os- 
cillations of the inter-exciton coherence in the ground-state of 
the vibrational mode \pxo-Yo{t)\ (thin red curve), c, Mandel 
Q-parameter of the vibrational mode for biological electronic 
coupling and d for zero electronic coupling. 



ten in the basis of exciton-vibration states of the form 
\X,n) becomes p(t ) = J2n p th{n) \X, n) {X, n\. Here 
n denotes the phonon occupation number of the collec- 
tive mode while Pth{n) denotes the thermal occupation 
of such level. We investigate the population of the lowest 
excitonic state /9yy(t) = 5^ n (y,7i|/o(t)|Y,7i), the abso- 
lute value of the coherence pxo-Yo{t) — (X,0\p(t)\Y,0) 
which denotes the inter-exciton coherence in the ground- 
state of the collective vibrational mode, and the non- 
classicality given by negative values of Q(t). Hamilto- 
nian evolution generates coherent transitions from states 
\X,n) to \Y,n + 1) (see Figure 1(a)) with a rate / 
that depends on the exciton delocalization (|V|/Ae), the 
coupling to the mode g, and the phonon occupation n 
i.e. / ~ g(2\V\/Ae)y/(n+l)/2. Since uj vih » K B T 
the ground state of the effective mode is largely popu- 
lated, such that the Hamiltonian evolution of the initial 
state is largely dominated by the evolution of the state 
\X, 0). This implies that the energetically close exciton- 
vibration state \Y, 1) becomes coherently populated at a 
rate / ~ g(2\V\/Ae), leading to the oscillatory pattern 
observed in the probability of occupation pyyit) as il- 
lustrated in Figure 1 (b). The low- frequency oscillations 
of the dynamics of pyy (t) cannot be assigned to the ex- 
citon or the vibrational degrees of freedom alone as ex- 
pected from quantum-coherent evolution of the exciton- 
plus-effective mode system. For instance, if the mode oc- 
cupation is restricted to at most n = 1, the period of the 
amplitude of pyy (t) is given approximately by the inverse 
of (l/2)(y/(AE - a) 2 + {2g 2 AV 2 /aAE)-(AE+a)) with 
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and (2g 2 AV 2 /a 2 AE 2 ) < 1. Coherent ex- 



citon population transfer is then accompanied by coher- 
ent dynamics of the inter-exciton coherence \pxo-Yo(t)\ 
with the main amplitude modulated by the same low- 
frequency oscillations of pyy (t) and a superimposed fast 
oscillatory component of frequency close to w V ib as can 
be seen in Figure 1 (b). This fast driving component 
arises from local oscillatory displacements: when V ~ 
the time evolution of each local mode is determined by 
the displacement operator with amplitude a(t) — 2g(l — 
exp(— ic<Mbi))/kMb • As the state \Y, l)(Y, 1| is coher- 
ently populated, the collective quantized mode is driven 
out of equilibrium towards a non-classical state in which 
selective occupation of the first vibrational level takes 
place. This manifests itself in sub-Poissonian phonon 
statistics as indicated by negative values of Q(t) shown 
in Figure 1(c). Similar phenomena have been described 
in the context of elec tron transport in a nanoelectrome- 
chanical systenPSMl. Importantly, such sub-Poissonian 
statistics arises only when the electronic interaction is fi- 
nite. For comparison, Figure 1(d) shows that if V = 0, 
a local excitation drives the mode towards a classical 
thermal displaced state with super-Poissonian statistics 
(Q(t) > 0). In short, non-classicality of the relative dis- 
placement mode quasiresonant with the excitonic tran- 
sition arises through the formation of exciton-vibration 
states. 



Energy and coherence dynamics - Interplay between 
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vibration-activated dynamics and thermal fluctuations 
leads to two distinct regimes of energy transport. For 
our consideration of weak electronic coupling, the co- 
herent transport regime is determined approximately by 
y/JTT c < 2g\V\/Ae. Here population of the low-lying cx- 
citon state is dominated by coherent transitions between 
exciton-vibration states and the rapid short-time growth 
of pyy{t) in this regime can be traced back to coherent 
transfer of population from \X, 0){X, 0| to \Y, 1}(Y, 1|. At 
longer time scales thermal fluctuations induce incoherent 
transitions from |X,0)(X,0| to \Y,0)(Y,0\ with a rate 
proportional to \/Afi C) thereby stabilizing population of 
\Y, 1)(Y, 1| to a particular value (see Figure 2(a)). In 
contrast, for \/\il c > 2g\V\/Ae population transfer to 
pyyit) is incoherent although the dynamics is still ac- 
companied by generation and coherent evolution of ex- 
citonic coherences (see Figures 2(g) and 2(h)). In this 
incoherent transport regime pyy{t) has a slow but con- 
tinuous rise reflecting the fact that the incoherent transi- 
tions from |X,0)(X,0| to |Y,0)(Y,0| now contribute sig- 
nificantly to the transport. However, transfer to pyy {t) 
is always more efficient with the quasiresonant mode than 
in the situations where only thermal-bath induced tran- 
sitions are considered (see dashed lines in Figures 2(d) 
and 2(g)). The underlying reason is that the system is 
transiently evolving towards a thermal configuration of 
exciton-vibration states. Hence, in both regimes excita- 
tion transfer to the lowest exciton state is achieved by 
transiently sharing energy with the relative displacement 
mode. The transition from coherent to incoherent ex- 
citon population dynamics is also marked by the onset 
of energy dissipation of the exciton-vibration system as 
shown in figures |2fc)j2ff) and ^i). While population 
transfer is coherent, energy dissipation into the thermal 
bath is prevented (Fig. ^c) and[2jf)). 

The beating patterns of the coherence pxo-Yo(t) re- 
veals the structured nature of the exciton-phonon inter- 
action and witnesses whether there is coherent exciton- 
vibration evolution. Because of this the frequency com- 
ponents of such oscillatory dynamics vary depending of 
the coupling to the thermal bath as can be seen in 
Figures [2jb), (e) and (h) . In the coherent regime 
Pxo-Yo(t) follow exciton populations with the main am- 
plitude modulated by the same relevant energy differ- 
ence between exciton-vibration states (see Figures J2jb) 
and (e)). In contrast, in the regime of incoherent trans- 
port the short-time oscillations of pxo-Yo(t) (between 
t = and t = O.lps) arise from electronic correlations 
due to bath-induced renormalization of the electronic 
HamiltoniarF^ and as such, the frequency of these short- 
time oscillations correspond to transitions of renormal- 
ized excitons. As the main amplitude is damped, this 
exciton coherence retains the superimposed driving at 
a frequency u; V ih>. The dynamical features presented in 
Figures 2(g) and 2(h) confirm previous findings based 
on perturbative calculations^ and agrees with the time- 
scales of the exciton coherence beating reported in the 
two-dimensional spectroscopy of cryptophyte algae^. 



Non-classicality - Interaction with a large thermal en- 
vironment leads to the emergence of classicality. In- 
deed, as the coupling of the electronic degrees of free- 
dom to the thermal background is increased, the maxi- 
mum non-classicality of the collective mode as indicated 
by the most negative value of Q(t) is reduced. How- 
ever, the time interval for which the mode remains in 
a non-classical state is not a monotonic function of A. 
For moderate values of A, the relative displacement mode 
spends longer periods in states with non-classical fluctu- 
ations i.e. periods for which Q(t) < as seen in Figures 
|3ja)-(b). Hence relatively medium coupling to a ther- 
mal background assists selective population of state \ Y, 1) 
thereby stabilizing non-classicality at a particular level. 
An important feature of the results presented in Figures 
3(a)-3(c) is therefore that non-classical features surviving 
the picosecond time scale are present across the whole 
range of thermal bath couplings. This is consistent with 
a slow decay of the exciton-vibration coherence pxo,Yi 
(not shown). Interestingly, the non-classical properties of 
the collective mode resemble non-classicality of thermal 
states (completely incoherent states) that are excited by 
a single quanta 41 . The non-classical properties predicted 
by Q(t) also agree with the non-classicality predicted by 
a parameter proposed by KlyshkxP^ as is described in 
appendix B. 

Non-classicality is witnessed by exciton population dy- 
namics. In order to illustrate this, we now investigate 
quantitative relations between non-classicality and en- 
ergy transport by considering relevant integrated aver- 
ages in the time scale of coherent exciton-vibration dy- 
namics denoted by r. For the parameters chosen, this 
time scale is about half a picosecond and corresponds 
to the time in which excitation energy would be dis- 
tributed away to other chromophorcs 16 . The time in- 
tegrated averages over r are defined as: (F[p(t)]) T = 
i jl dt F[p(t)] , where F[p(t)] corresponds to the exciton 
population pyy (t) and the non-classicality of the collec- 
tive quantized mode through periods of sub-Poissonian 
statistics Q(t)Q[— Q(t)] as functions of the coupling to 
the bath A. As shown in Figure [4] the average exci- 
ton population and non-classicality follow a similar non- 
monotonic trend as a function of the coupling to the 
thermal bath, indicating a direct quantitative relation 
between efficient energy transfer in the time scale r and 
the degree of non-classicality. The appearance of max- 
imal point in the average non-classicality as a function 
of the system-bath coupling indicates that the average 
quantum response of the relative displacement mode to 
the impulsive electronic excitation, is optimal for a small 
amount of thermal noise. 

The above results indicate that the degree of dereal- 
ization of the initial exciton state is key to enabling non- 
classical fluctuations of the collective quantized vibra- 
tion. This suggests that statistical mixtures of excitons 
with finite purity can still trigger such non-classical re- 
sponse. To illustrate this we now consider mixed ini- 
tial states of the form p(to) = g cyi ® g^ h where g cyi — 
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FIG. 2. The dynamics of pyy(t) (top row), |pxo,yo(£)| (middle row) and Energy of the exciton vibration system (bottom row) 
for three interaction strengths to the low-energy thermal bath: a-c: A = 6cm" 1 , d-f A = 20cm -1 and g-i: A = 110cm -1 
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FIG. 3. a-c Dynamics of the Mandel Q-parameter for A = 6, 20, 110cm 



r \X) (X\ - (1 - r) \Y) (Y\ with 1/2 < r < 1. The as- 
sociated linear entropy quantifying the mixedness of the 
initial exciton states is given by Sl — 2r(r — 1). The 
time-averaged non-classicality (Figure [5] (a)) and aver- 
age population transfer (Figure [5] (b)) follow similar de- 
creasing monotonic trends for mixed states. These find- 
ings then suggest direct quantitative relationships be- 
tween non-classicality and exciton population transport 
in the relevant time scale. Outlook.- Our work shows 
that in molecular dimers operating in the weak electronic 
coupling regime, exciton-vibrational resonances can trig- 
ger a non-classical response of a quasiresonant collective 
high-energy vibration thereby taking advantage of quan- 



tum fluctuations to efficiently populate a target exciton 
state. These non-classical fluctuations are observed for 
a variety of initial exciton conditions including statisti- 
cal mixtures suggesting that such non-classicality may 
be triggered even under incoherent input of electronic 
population. These quantum phenomena could be ex- 
perimentally accessible by transient coherent ultrafast 
phonon spectroscopy^ which is sensitive to low-phonon 
populations of high-energy vibrations^. The framework 
here discussed could also be used to give more insight 
into the possible non-classical response of vibrational 
modes in other chemical sensors^ which are conjec- 
tured to operate through weak electronic interactions^ 
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mechanism for thermodynamic control. Understand- 
ing how quantum phenomena affect the thermodynamics 
of energy conversion in biomolecules could provide an 
unifying framework to quantify the advantages of quan- 
tum phe nome na in a variety of biomolecular transport 
processe d 48 1 49 1 . 
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FIG. 4. Time integrated averages of exciton popu- 

lation pYy(t), inter-exciton coherent in the lowest vibra- 
tional level \pxo-Yo{t)\ and non-classicality as quantified by 
Q(t)6[—Q(t)] (blue, red and green respectively) as functions 
of coupling to the thermal background by fixing environment 
cut-off frequency S7 C = 100cm -1 and varying reorganization 
energy A. 
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FIG. 5. The time integrated averages of exciton popula- 
tion pYy(t) — pyy(0) and non-classicality as quantified by 
Q(t)Q[—Q(t)] as functions of the linear entropy Sl = 2r(r — 1) 
for < r < 1/2. Curves are for different system-bath cou- 
plings. 



to sense vibrational frequencies of the order of a thousand 
wavenumbers^. 

From a more general perspective, our work illustrates 
how exciton-vibrational resonances allow for rich dynam- 
ics of energy and coherences thereby providing a mech- 
anism for adjusting the rates of entropy production and 
output energy. In this context, the associated vibrational 
dynamics could be understood as an internal quantum 
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Appendix A: Hierarchical expansion of exciton-vibration 
dynamics 

Here we provide details of the model for a system 
consisting of N electronic degrees of freedom each in- 
teracting with a quantized vibrational mode and the 
method used to accurately calculate the dynamics of 
these exciton-vibration systems. The system Hamilto- 
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nian Hg = H c \ ® l vib + l c i ® H vih + H c i_ vih , is defined 
by 

N N 

h c1 = E ^t°T + E v * + 4°t ) ( A1 ) 

i—X * , j < i 



H, 



vib 



i=l 



and 



H, 



ol-vib — gi&t a 



(A2) 



(A3) 



Here, the operators af 1 create/annihilate electronic exci- 
tations at site i, while b\ (bi) create (annihilate) phonons 
associated with the vibrational mode of site i. Si denotes 
the excitation energies and the strength of interac- 
tions between sites. u;^ ib is the frequency of the vibra- 
tional mode at site i and gi = -\/5Wvib the coupling of 
the mode to its site. Si is the Huang-Rhys factor or 
mean number of phonons in the 'polaron cloud' formed 
when the mode is fully displaced. The the high energy 
modes considered here, we expand the quantized vibra- 
tional modes of the system in the basis of single phonon 

(Fock) states (e.g. b t = £f (n + l ) 1/2 \ n ) ( n + M) trun- 
cating at M phonons. This system is then coupled to a 
continuum of bosonic modes (the thermal background) 



Hb = Ek w k^k, with an interaction 



H 



i = EE 9k,i(trt"i ® 1 vib)(&i c + b k 



(A4) 



such that a mode with frequency Wk is linearly coupled 
to the electronic excitation at site i with a strength gk,i 
as specified by a spectral density Ji(co) — £ k \gk,i\ 2 5(u} — 
Wk). We assume the bath is spatially uncorrelated and 
identical for each site. 

An infinite hierarchy of coupled differential 



equation; 



is used to express the dynamics of 



the exciton-vibration 
These read as, 



system density matrix ps(f). 



N K 



d t p n (t) = -i/h[H s ,p n {t)] -^2^2n lk v ik p n (t) 

i=l fc=0 
N / oo K \ 

t \ , j ^ik , t ^ik I 

4=1 \fc=l k=l / 

JV K 

~ *EE V ( n ik + l)|Cjfc| 



i=l k=0 
N K 



for the case of a Drude spectral density and in a rescaled 
formulation^ which admits their efficient adaptive time- 
step numerical integration. Operators p n are indexed by 
the multi-index n which has entries no, and n — n%k- 



n ik 1S n w ith entry n 



jk 



± 1. The reduced den- 



sity matrix of the system p$ is given by the operator 
po. The coefficients Vi k and Ci k are the Matsubara fre- 
quencies and coefficients appearing in the (truncated) ex- 
ponential decomposition £fc=o c ik&~ Uikt of the bath cor- 
relation function applied in ReP^. In the above terms 
k > K are truncated with a Markovian approximation 
scheme^. This results in accurate dynamics for K = 1 
at high temperatures (/3f2 c < 1). In the present exciton- 
vibration case the system-bath coupling operator takes 
the form Qi = afa~ ® l v ib- 



Appendix B: B parameter 



Modulations of adjacent phonon number occupation 
as quantified by negative B n , guarantees negative re- 
gions of a quasi-probability distribution^ as is a poten- 
tially an experimentally accessible quantity. The non- 
classicality indicated by sub-Poissonian statistics is con- 
firmed by negative B = 2P(0)P(2) - P(l). This arises 
from selective population of the P(l) level of the relative 
displacement vibrational mode. 
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